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We investigate the conformations of DNA-like stiff chains, characterized by contour length (. L ) 
and persistence length (Z p ), in a variety of crowded environments containing monodisperse soft 
spherical (SS) and spherocylindrical (SC) particles, mixture of SS and SC, and a milieu mimicking 
the composition of proteins in E. Coli. cytoplasm. The stiff chain, whose size modestly increases 
in SS crowders up to <f> « 0.1, is considerably more compact at low volume fractions (0 < 0.2) in 
monodisperse SC particles than in a medium containing SS particles. A 1:1 mixture of SS and SC 
crowders induces greater chain compaction than the pure SS or SC crowders at the same <f> with 
the effect being highly non-additive. We also discover a counter-intuitive result that polydisperse 
crowding environment, mimicking the composition of a cell lysate, swells the DNA-like polymer, 
which is in stark contrast to the size reduction of flexible polymer in the same milieu. Trapping 
of the stiff chain in a fluctuating tube-like environment created by large-sized crowders explains 
the dramatic increase in size and persistence length of the stiff chain. In the polydisperse medium, 
mimicking the cellular environment, the size of the DNA (or related RNA) is determined by L/l p . 
At low L/lp the size of the polymer is unaffected whereas there is a dramatic swelling at intermediate 
value of L/lp. We use these results to provide insights into recent experiments on crowding effects 
on RNA, and also make testable predictions. 

The recognition that the crowded cellular environ¬ 
ment can profoundly influence all biological processes, 
such as gene expression m, protein HI] and RNA 
folding mm and protein-protein interactions LU • is 
getting increasing attention recently although its im¬ 
portance was recognized decades ago m- A simple 
calculation using the typical concentration of macro¬ 
molecules shows that the average spacing between pro¬ 
teins in the E. coli is ^ 4 nm, which is comparable to 
the size (radius of gyration) R g « 0.37V 1//3 nm of a 
folded protein m with N ~ 300 amino acid residues. 

Therefore, the cellular interior, replete with macro¬ 
molecules of different sizes and shapes (a polydisperse 
soup), is crowded, affecting the stability and shapes 
of the molecules of life. For example, compaction 
of DNA, relevant in a variety of biological processes 
ranging from organization of nucleoid in bacteria to 
DNA packaging in phage heads, is greatly facilitated 
in the presence of neutral osmotic agents. The effect of 
neutral polymer (poly-ethylene oxide) (PEO) in com¬ 
pacting DNA (T-condensation) was demonstrated in 
a pioneering study by Lerman PH. who showed that 
DNA undergoes a dramatic reduction in size if the 
concentration of PEO exceeds a critical value. Be¬ 
cause the interactions between DNA and PEO in these 
experiments was established to be repulsive [13], 4/- 
condensation is determined solely by the volume ex¬ 
cluded to DNA by the crowding particles. Subse¬ 
quently, Post and Zimm m produced insightful the¬ 
oretical explanations based on further experiments on 
crowding-induced compaction of DNA. In the inter¬ 


vening years a number of theoretical and experimental 
studies have explored various aspects of DNA com¬ 
paction SMI]. 

Despite these advances a molecular description of 
how macromolecular crowding, especially a milieu con¬ 
taining polydisperse crowders, affects the spatial orga¬ 
nization of DNA is poorly understood. Simulations 
of DNA in the presence of explicit crowding particles 
are computationally intensive, but their importance 
in describing the shapes of flexible polymers has been 
demonstrated in a number of recent studies p~8l [22] . 
Because of technical complications theories do not take 
into account the effects of polydispersity mnn [13- 
i30q In addition, the extent of structural changes and 
how they depend on mixtures of crowding particles of 
different shapes for realistic sizes of crowders found in 
the cytoplasm are unknown. Most of the studies have 
focused on intrinsically flexible polymers with very lit¬ 
tle work on relatively short stiff DNA-like polymers, 
and related RNA. Stiff short chains can display quan¬ 
titatively different behavior than the one associated 
long flexible polymers, as argued in the case of cy- 
clization of DNA [3TH33] . These issues take on added 
importance because of potential relevance to genome 
confinement |4] 34, 35 , where the structural organiza¬ 
tion is functionally related to gene expression. 

Inspired by these observations, we first performed 
simulations of worm-like chain [36] (WLC - a reliable 
polymer model for describing many of the properties 
of DNA at high salt concentration) in monodisperse 
solutions containing soft spherical (SS) and sphero- 
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cylindrical (SC) crowding particles. Here, we focus on 
DNA-like chains in the limit where the contour length 
(L = 7V m / 0 , N m is the number of monomers, and Iq 
is the bond length) is not significantly longer than the 
persistence length, l p . The systematic study leads to a 
number of unexpected predictions, which can be tested 
using synthetic polymers in the presence of nanoparti¬ 
cles and DNA using crowding agents, (i) The SC Crow¬ 
ders induce greater compaction than SS particles. The 
compaction is accompanied by substantial reduction of 
l p (by nearly a factor of two) similar in magnitude to 
that observed by ion-induced compaction of DNA and 
RNA [37] . (ii) One of the most striking results of our 
study is that stiff chains are more compact in a mix¬ 
ture of SS and SC, at physiologically relevant volume 
fractions, than in monodisperse crowders at the same 
volume fraction. The substantial compaction of DNA 
in the mixture is due to ordering (dense packing) of the 
SS crowders around DNA due to depletion attraction 
induced by the SCs. (iii) We also carried out sim¬ 
ulations of DNA in polydisperse crowders mimicking 
the composition of macromolecules in E. coli. Sur¬ 
prisingly, we obtained a counter-intuitive and theoret¬ 
ically unanticipated result that polydisperse spherical 
crowding particles cause swelling of DNA at volume 
fraction « 0.3 appropriate for E. Coli. cytoplasm. Sur¬ 
rounded predominantly by large-sized crowders, the 
swollen conformation of DNA (nearly three times in¬ 
crease in volume of the chain) is entropically driven. 


RESULTS 

Crowding-induced softening of DNA. In order 
to set the scale for bending energy of the polymer, we 
first calculated the end-to-end distance distribution, 
P(R ee ), in the absence of crowders to obtain the 
bare persistence length, l p for (j) = 0. By fitting the 
simulated P (r) to an analytic expression for WLC 
(Eq.S3) [3L, we obtained l p = (15.4 ± 0.08) a m ~ 49 
nm using cr m ~ 3.18 nm (fib), which coincides with 
l p for DNA in high monovalent salt concentration. 
Because / p /cr m > 1, the shape of the chain in the pres¬ 
ence of crowding particles should be determined by an 
interplay of bending rigidity and attractive depletion 
interaction due to crowders. P(R ee ) of the WLC with 
varying for both monodisperse SS and SC crowders 
show gradual shift to the smaller R ee with increasing 
(j) from 0 to 0.4 At high </>, as 0 increases, l p 

obtained from the fit of P(r) to Eq.S3 gradually 
decreases for both SS and SC crowders, implying 
that the crowding particles induces compaction of 
the WLC polymer ([TJd, top panel). Interestingly, the 
WLC polymer exhibits a non-monotonic dependence 



FIG. 1: (a) Coarse-grained models of spherical (SS) Crow¬ 

der (red), spherocylinderical (SC) crowder (yellow) and 
WLC chain (blue) used to model DNA. The relevant di¬ 
mensions are labeled <7 S phi cr C yi , and cr m . (b) End-to-end 
distance distribution of WLC at two different volume frac¬ 
tions (0 = 0.2, 0.4) of SS and SC crowders. The top panel 
shows the persistence length (l p ) of the polymer as a func¬ 
tion of the volume fraction of SS and SC crowders. 


of size with increasing 0. The SS crowders induce a 
minor increase of l p (stiffening or expansion) of the 
chain & top panel and|2^i) for 0 < <p < 0.1, followed 
by a decrease of l p (softening or compaction) at larger 
4> = 0.2 — 0.4. In contrast, the SC crowders reduces 
l p much more efficiently than the SS crowders at 
< j) < 0.3 in that an expansion similar to the one in the 
SS crowders at 0 ^ 0 — 0.1 is not observed. Instead, 
there is a modest re-stiffening of the DNA due to 
SC crowders when </> is in the range from = 0.2 to 
< j) = 0.4. We further substantiate this result below by 
calculating the change in the polymer size (R g ) and 
nematic order parameter of the crowders. The slight 
increase in l p not withstanding, the overall trend is 
that there is substantial softening (, l p decreases by 
nearly a factor of 2) as 0 increases from 0 to 0.4. 

Dependence of R g on 0 for monodisperse 
crowders. Snapshots of polymer conformations at 
different values of </> show modest compaction as (j) in¬ 
creases (|2r). Similar to Z p , the dependence of R g on </> 
for SS and SC crowders displays substantial difference. 
At = 0.2, R g { 0.2) is smaller than R g ( 0) by only 4 
% in SS crowders whereas R g ( 0.2) decreases by 17 % 
in SC crowders. Given that the volume of the chain 
is ~ R g {(t)) the extent of compaction induced by SC 
is substantial compared with SS crowders. The quan¬ 
titative difference between the effects of SS and SC 
crowders on the chain compaction is explained using 
the depletion interaction (or Asakura-Oosawa (AO) in¬ 
teraction [38] ) that produces an effective attraction be¬ 
tween monomers. The strength of the AO interaction 
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FIG. 2: (a) The change of R g of the polymer from R g (0) = 37.6 nm with increasing 0. (b) Liquid crystal order parameter 

for SC crowders, averaged over the ensemble of crowder particles, as a function of distance, (c) The extent of local nematic 
ordering is quantified with Q 2 (r) at r — 4.2cr m . 


for SS, roughly given by ~ (fiksT/cr^ ph , has to exceed 
the energy ^(P/Z p )&bT) required to bend the poly¬ 
mer on scale l p for compaction to occur. For small 0, it 
is unlikely that the AO attraction can compensate for 
the bending penalty. Thus, we expect little change in 
R g ( (j )) at small </> for SS crowders. On the other hand, 
the strength of the AO interaction on the WLC for 
SC crowders is ~ (f)P(a rn /al yl )kBT where P(= 2 <J cy i) 
is the cylinder length. In both cases, the origin of 
the AO depletion interaction, which has to exceed the 
bending energy to compact the stiff chain, leading to 
an effective short range (on length scale ^ cr m ) attrac¬ 
tion between monomers and polymer compaction, is 
purely entropic. [2 ^l shows that R g (<fi) decreases mono- 
tonically till 0 ~ 0.3 for SC, indicating that for the 
parameter used, the strength of the attractive AO in¬ 
teraction due to SC crowders can exceed the penalty 
for bending the chain on scale ^ l v . 

The strong effect of compaction of the WLC chain 
induced by SC relative to SS crowders can also be 
quantified by comparing the volume excluded to the 
polymer by the crowders. On the basis of scaled parti¬ 
cle theory, we can estimate the entropy cost of insert¬ 
ing a hard sphere of dimension cfhs in a box containing 
hard fluid particles. The entropy difference for insert¬ 
ing the hard sphere of diameter cfhs is related to C(a) 

ED, 

x _ Kyi _ 3 [(cfHS + &cy\) 2 P + {&HS + ^cyl) 3 ] 

Py&HS) 77 ~ “3 

Lsph 2 [CFHS + cr S ph) 

(i) 

where V cy i ( V sp h ) is the volume excluded by rod-like 
(spherical) crowders. For the parameters listed in Ta¬ 
ble SI, we find that V cy i > V sp h provided cfhs ~ R g > 
a cy i or cr sp h- The entropic cost of inserting a spherical 
particle of size R g into a fluid of cylindrical crowders 
exceeds that for inserting it into a system consisting of 
spherical crowders. By achieving greater compaction 


of the WLC in SC crowders, the entropy difference is 
minimized, thus, explaining the results in lit- 

The pattern of compaction in the two crowd¬ 
ing environment is qualitatively different. For SS 
crowders, compaction occurs only when 0 exceeds 
~ 0.1 ( dR g /d(j) « 0 as (j> —>> 0). Although the trend 
is not so clear as in l p ((f)) top panel), R g ((l)) 

&) shows a signature of minor swelling in the 
range of </> = 0 — 0.1, which was also observed in 
experiments on the effects of crowding on a ribozyme 
(see below) [8]. In contrast, for anisotropic crowding 
agents (SC crowders), R g (0) decreases monotonically 
for (p < 0-3 ( dRg/dcj) < 0 for all </>) and increases 
from <f> = 0.3 to 0 = 0.4. High volume fraction 
of SC results in the reswelling of WLC (Figure 
2a), showing that the shape of the crowding par¬ 
ticles can have a profound effect on size of a stiff chain. 

Local nematic order and increase in R g . In¬ 
terestingly, as (p exceeds 0.3, R g {4>) of the chain in¬ 
creases in SC crowders §t). The increase in R g ((j)) 
at higher </> is due to plausible development of lo¬ 
cal nematic ordering for SC crowders at <p > 0.3. 
The system consisting of pure SC crowders will un¬ 
dergo an isotropic to nematic phase transition if (j) ex¬ 
ceeds a critical value We calculated the liq¬ 

uid crystal order parameter, (P 2 (cos#)), where 0 is 
the angle between the long axes of any pairs of SC 
crowders, P 2 (x) = \ (3x 2 — l) is the second Legen¬ 
dre polynomial, and (...) denotes the ensemble av¬ 
erage. (P 2 (cos $)) is almost zero (< 0.05) for all </>, 
which means that even the highest <p(= 0.4) is less than 
However, locally the crowders adopt nematic¬ 
like state. To ascertain if this is the case, we calcu¬ 
lated P 2 (r) = (^2ijP2 (cos 6ij)S(\fi -fj\ -r)) as a 
function of distance, r. It shows that the angular 
correlation of SCs becomes stronger as (j) increases at 
short distances. 
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To quantify the extent of “local” nematic ordering, 
we calculated Q 2 ( 7 *) = Np 1 f r P 2 (r / )dr / , where 
Np is the number of pairs separated by r and r m i n 
is the minimal distance between SC crowders. At 
the distance r = 2.3 a m where the pair correlation 
has first peak for 0 = 0.4 it), P2(2.3cr m ) increases 
significantly from 0.0 to 0.23 as 0 increases from 0.1 
to 0.4 if)- At r = 4.2cr m , Q 2 W ~ 0.23 suggests that 
the local nematic ordering of cylindrical crowders is 
reinforced in the vicinity of the WLC. This result 
implies that relatively stiff chain induces ordering of 
SC crowders along the polymer axis, and strengthens 
the anisotropic interaction of rod-like particles. At 
high 0 the chain may be thought of as being in a local 
nematic field, which elongates the polymer along the 
local direction of the nematic field, thus explaining 
the increase in R g {4>) when 0 exceeds 0.3. 




FIG. 4: (a) RDFs are calculated for all possible combi¬ 

nations of crowder-crowder pairs (at (i) 0 = 0.2 and (ii) 
0.4) and (iii) monomer-crowder (at 0 = 0.4). Distribution 
of SS and SC crowders around the WLC chain is depicted 
using a snapshot from the simulations (The particle sizes 
are not in scale. Although cr sp h — 4cr m and a cy i = 4 2 ^ 3 cr m 
were used in actual simulations, we deliberately reduced 
the sizes of SS and SC crowders for clear illustration of 
the crowders around the chain), (b) RDFs of SS-SS and 
SC-SC pairs at 0 = 0.1 — 0.4 in monodisperse condition 
are presented to highlight the enhanced local ordering of 
SS and SC crowders in the mixture. 


FIG. 3: Compaction of stiff chain in a mixture of SS and 
SC crowders. R g of WLC as a function of 0 of 1:1 SS and 
SC mixture (green). The orange line is the calculated R g 
at each 0 value by assuming that the effects of SS and SC 
crowders on WLC compaction are additive. The results of 
monodisperse SS and SC crowders in[2^ are shown with the 
dashed lines to underscore the substantial enhancement of 
the chain compaction by the mixture. 

Non-additive effect in a mixture of spheres 
and spherocylinders. The dependence of R g (0) on 
0 of the WLC in the 1:1 mixture of SS and SC crowders 
is shown in[3j The mixture has a profound effect on the 
size of DNA. The value of R g (0.4) is reduced by over 
40% from R g ( 0), whereas the maximum compaction in 
monodisperse SC at 0 = 0.4 is only 10%. In order to 
highlight the non-additive effects of the two crowders 
with different shape, we also show, at 0 = 0.2 and 0 = 
0.4, the expected result for i?^(0) = R g ( 0) + SRf (0), 
where AR^(0) = AR^ s (0/2) + 8R? g c {(j)/ 2), obtained 
using the assumption that the effect of crowders on the 
extent of R g reduction is additive. Remarkably, R g (4>) 
is significantlly lower than R^{4>) ([ 3 ]), (or \SR g ((j))\ > 
|AR^(0)|, |ARS c (0)|, |AR^ s ( 0)|), indicating that the 


mixture of SS and SC restricts the volume available to 
the WLC to a much greater extent than the individual 
components do. 

The surprising finding of significant compaction in 
the mixture can be qualitatively explained using the 
notion of depletion potentials for mixture of SS and 
SC. Consider the interaction between two spherical 
particles in the presence of rods (SC particles). If a SC 
particle is spatially trapped then the rod loses transla¬ 
tional and rotational entropy because of orientational 
restrictions. The large unfavorable entropy loss results 
in the depletion of the SC from the space where the SC 
particle is trapped. The result is that there would be 
an excess osmotic pressure due to the AO attraction 
that pushes the SS and the monomers together. Two 
consequences of the entropy-driven depletion interac¬ 
tions are : (i) Due to the attractive interactions, the SS 
particles are more closely packed than in the absence of 
the SC (compare #™ss( r ) a t 0 from[4^ and #ss°ss( r ) 
at 0/2 from & (monodisperse crowders)). Indeed, in 
the limit of a cy i < P < g sv h^ it has been experimen¬ 
tally shown that the addition of a small fraction (by 
volume) of the SC can even lead to crystallization of 
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low density suspension of hard spheres mm • (ii) 
We also expect that the excess volume available to 
the WLC should be greatly reduced compared to the 
monodisperse crowders. In such a confined space the 
WLC should be considerably more compact than in 
the presence of monodisperse crowders at the same <fi. 

The expected enhancement in the packing of the SS 
due to the SC is evident in the pair distribution func¬ 
tion g(r). The results for 0 = 0.2 and cj> = 0.4 
(i), (ii)) show that both the radial distribution func¬ 
tions (RDFs), <7sc* sc (r), between SCs and g^sc ( r )> 
between SC and SS, do not exhibit significant struc¬ 
ture. In sharp contrast, g^s^ss ( r ) has the structure 
corresponding to a high density liquid (especially at 
0 = 0.4 in[4^-(ii)), which is remarkable given that at 
4>/2{= (f>ss = 0-2), corresponding to the same volume 
fraction occupied by SS, g™s~ss( r ) (jU 3 ) * s relatively fea¬ 
tureless. 

The much stronger depletion force due to the 
SC crowders results in a considerable reduction in 
the volume accessible to the chain, which explains 
the dramatic reduction of R g (</>) compared to the 
pure component case. The WLC at = 0.4 is 
surrounded predominantly by SS crowders, which is 
shown by g™!x( r )> tlie RDF between the monomer 
and the crowders (X = SS or SC) |4^-(iii)). The 
number of SS near the monomer calculated using 
N S s = 47r (^) / Q rmi " r 2 g m _sc (r) dr where r min is the 
first minimum in the g m -ss (0 at = 0.4 is 11.4. A 
similar calculation for SC gives ATsc = 4.7. 

Polydisperse crowders mimicking the E. coli 
cytoplasm composition. To a first approximation 
the E. coli cytoplasm may be represented by a mixture 
of spheres because majority of the crowders present in 
large numbers (ribosomes, polymerases and other large 
complexes, and smaller particles) are compact [42] . In 
order to assess the shape of DNA-like chain in such 
a mixture, we investigated the effect of polydisperse 
spherical particles on the conformational fluctuations 
of the stiff polymer. 

Strikingly, the behavior of the stiff chain in a poly¬ 
disperse mixture of SS particles differs drastically from 
those in the mixture of SS and SC of the same volume 
<i- A few features in the non-monotonic dependence 
of R g (0) as a function of 0 &) are worth pointing 
out. (i) There is a very modest reduction in R g (0) at 
0 ^ 0.1, which has only large crowding particles. Such 
a decrease is comparable to that found in[2^. (ii) Unex¬ 
pectedly, R g ((/)) starts increasing in a mixture contain¬ 
ing large and medium sized crowders. Most strikingly, 
in the mixture roughly mimicking the E. coli cyto¬ 
plasm [42] there is a large increase (^ 40%) in R g (</>) 
compared to 0 = 0. The competition between bending 


free energy and depletion potential leading to a dra¬ 
matic swelling of the stiff chain is counter-intuitive. 
The ensemble of the chain conformations (j5^i) exhibit¬ 
ing the expansion of the chain captures these effects 
visually, (iii) The E. coli mixture dramatically stiffens 
the polymer. The persistence length of the chain for 
the swollen chain in the E. coli mixture is about 2.6 
times larger than for the one at <fi = 0. The stiffening 
effect of mixture of spherical crowders on WLC cap¬ 
tured the snapshots in [5] (iv) It is noteworthy that 
in the polydisperse mixture with E. coli. composition 
reduces the size of a flexible chain that lacks bending 
penalty & and Figure S3), underscoring the impor¬ 
tance of chain stiffness. 

The reswelling at high (f) &#) can be understood 
qualitatively using the following arguments. At a spec¬ 
ified total volume fraction there are a lot more small 
crowding particles than large ones. Therefore, the en¬ 
tropy of the crowding particles is maximized if the 
WLC is surrounded by the larger sized particles with 
the smaller ones being further away from the chain. In 
this picture, the WLC is localized in a region in which 
the larger sized particles are with higher probability in 
proximity to the monomers (Figure S4). Because the 
interactions between the crowding particles and the 
stiff chain are repulsive the DNA chain would prefer to 
be localized in a largely crowder-free environment. If 
we assume that such a region is roughly spherical, cre¬ 
ated predominantly by the largest crowders, then its 
size has to be on the order of R? g to accommodate the 
WLC chain. In such a cavity there is an entropic cost 
to confine the stiff chain. The probability of finding 
such a region decreases exponentially as R g gets large. 
In addition, in a spherical region the DNA would form 
spool-like structures requiring overcoming bending en¬ 
ergy. The combination of these effects makes it likely 
that an optimal spherical regime can be found to min¬ 
imize the free energy of WLC. If the region is cylin¬ 
drical and large enough such that tight hairpins (cost¬ 
ing substantial bending penalty) are avoided then the 
chain free energy may be minimized by confining it in 
a roughly cylindrical cavity. Such a possibility is sup¬ 
ported by simulations, which show that on an average 
the shape of the depletion zone is aspherical resembling 
a fluctuating tube (see images in [6^i) . As a result, we 
can visualize the polymer to be essentially confined to 
an anisotropic (but fluctuating) tube in which trans¬ 
verse fluctuations of the chain are restricted but one 
in which tight hairpin turns cannot form because of 
bending penalty. In such a cavity the chain stiffens, 
thus expanding in size in order to minimize both the 
bending penalty and entropy cost of confinement. A 
more quantitative and accurate theory is difficult to 
construct because of the many body correlation among 
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the poly disperse crowding particles. 

The plausibility of the physical picture given above 
is further substantiated by examining how the crowd¬ 
ing particles with different sizes are arranged in space 
and how they surround the extended WLC. Distribu¬ 
tions of the polydisperse crowding particles are not 
uniform, but exhibit local size ordering. This is evi¬ 
dent in the local size correlation function (r-dependent 
size variance), £(r) = (didj) r — (d) 2 where (d) is 
the mean diameter and (didj) r in the first term de¬ 
notes an average of the product of two diameters di 
and dj taken over crowders (j) located at a distance r 
from a crowder i [43 . The local size correlation £(r) 
(we set (didj} r = 0 when no pair exists), which is 
by definition zero for both monodisperse crowders and 
homogeneously distributed polydisperse crowders, re¬ 
veals the presence of size ordering up to r < 15 — 20 
nm. This implies that there is a non-uniform ordering 
in the mixture of spherical particles, which is respon¬ 
sible for reswelling ([5^). 

The fraction of volume occupied by three different 
crowders (</>x lx (r) ( = if ( 7^) 3 x 0 ^ 5 , m-x( r )) with 
X = 1,2,3) as a function of distance from the WLC 
monomers in |§ further captures the non-uniform dis¬ 
tribution of crowders. The crowders with intermediate 
size (X = 2 ) occupy the largest volume near the WLC. 
In addition, the comparision of </>™ lx (r) with 0 ™ ono ( r ) 
(the volume fraction of monodisperse 11 %-crowders 
around WLC monomer; the dashed line in [ 6 ^) shows 
that the intermediate (X = 2) and large (X = 1) sized 
crowders are pushed closer to the monomers by the 
small crowders, which confines a segment of the poly¬ 
mer to a tube-like region §6^). The depletion forces 
in a polydisperse solution give rise to a spatial inho¬ 
mogeneity of crowders, resulting in the chain being 
confined to a cylindrical region created by the large¬ 
sized (X = 1,2) crowders. The expansion of the chain 
in such a confined space [44] provides a plausible phys¬ 
ical explanation for the large increase in the size of the 
DNA. 


DISCUSSION 

Swelling and collapse of DNA. The counter¬ 
intuitive finding that a stiff polymer, with 7 =ss L/l p not 
large (see below), can swell relative in a polydisperse 
mixture of spheres, is (to our knowledge) unprece¬ 
dented. There are simulation and theoretical studies 
predicting the collapse of flexible polymers and pro¬ 
teins in mixed solvents due to volume exclusion effects 
alone [45l [46] . However, the present study shows pre¬ 
cisely the opposite behavior for stiff chains whereas a 
flexible chain tends to become compact (not a glob¬ 


ule in the E. Coli. like milieu). We propose that this 
unusual effect is related to an interplay of chain bend¬ 
ing and the complex depletion effect in a polydisperse 
crowding system. In order to substantiate our proposal 
we carried out simulations for chains with 7 varying 
from 2 to 10 in the model E. Coli -like system. The 
simulations show hardly any change in R g (see Figure 
S5). For all the values of y(< 10) the stiff polymer can 
be accommodated in a large enough region in which 
the crowding particles do not suppress the conforma¬ 
tional fluctuations. Only when 7 exceeds a minimum 
value, but is not too large, then chain swelling occurs 
by formation of tight turns. 

We provide arguments that when 7 exceeds a certain 
value the DNA-like polymer must undergo a transition 
from the swollen state to a globule. In other words, 
there must be a sharp coil-globule transition induced 
by the crowders. When 7 1 chain stiffness is not 

that relevant and the polymer behaves like a flexible 
polymer. In this limit, using our results in ref. m we 
predict that 7 has to exceed 40 to observe a genuine 
coil-globule transition. The simulation results and the 
physical arguments allows us to predict a rich depen¬ 
dence of R g in the E. Coli. environment @o- 
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FIG. 5: Effects of polydisperse soft sphere mixture on the 
size of WLC {L/l p — 20 at </> = 0). (a) R g is calculated 
for WLC (black solid line) and SAW (orange dashed line) 
polymers at (i) </> = 0.11 of crowders with n = 10.4 nm, (ii) 
(p = 0.22 with 7*1 = 10.4 nm (11 %) and 7*2 = 5.2 nm (11 %), 
and (iii) 0 = 0.3 with 7*1 — 10.4 nm (11 %), r 2 = 5.2 nm (11 
%), 7*3 = 2.6 nm (8 %). The largest crowder represents the 
ribosomes, the 5.2 nm crowders correspond to polymerases 
and other large protein complexes, and the smallest sized 
particles is the average size of other crowders in the milieu. 
Structural ensembles of polymers at </> = 0 and 0.3 are 
shown at the bottom, demonstrating the contrasting effect 
of polydisperse crowders on the conformations of WLC and 
SAW. (b) P(R g ) (top) and P{R ee ) (bottom) of WLC at 
4> = 0 and 0.3. The fit of P{R ee ) at (p = 0.3 to Eq. S3 
in the text gives l$ =0 ' 3 — 41 cr m ~ 130 nm, which is ~ 2.6 
fold greater than /p =0 (~ 15cr m 497i77i). P(R g ) of SAW 
(flexible self-avoiding chain) is shown in Figure S3. 
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FIG. 6: (a) A snapshot from simulation, demonstrating 

(left) the polydisperse crowding environment, (middle) the 
WLC inside the crowders, and (right) the crowding par¬ 
ticles decorating the monomers, (b) Local size ordering 
correlation function, £(r), indicates a non-uniform, hetero¬ 
geneous size ordering of polydisperse crowders. (c) 0™ lx ( r ) 
(X = 1, 2, 3), the fraction of volume around a monomer oc¬ 
cupied by three different types of crowders in polydisperse 
crowding environment (case (iii) in (a)) as a function of 
monomer-crowder distance (solid lines). The dashed line 
shows the corresponding quantity for monodisperse crowd¬ 
ing environment. 


Insights into crowding effects on RNA. Re¬ 
cent experiments have examined the effects of poly- 
tehylene glycol (PEG) on the folding of a ribozyme 
[8]. It has been argued that the impact of PEG can 
be understood based on the excluded volume effect. 
The SAXS experiments on Azoarcus ribozyme with 
195 nucleotides shows that R g , near the midpoint of 
the Mg 2+ ion needed for the folding transition, ini¬ 
tially increases before becoming compact ([ 7 ) 3 ) . Fold¬ 
ing of this RNA is also accompanied by decrease in the 
persistence length, which can be modulated by crow¬ 
ders. For Azoarcus ribozyme 7 ~ 36 where we have 
used L « 195 x 0.55 nm and l p ~ 3 nm [47]. If the 
theoretical prediction for the monodisperse SS in Fig¬ 
ure 2a is correct then we expect a modest increase 
in R g as PEG (assumed to be sphereical) concentra¬ 
tions increase. The experimental data is in qualitative 
agreement with this expectation. It would be most in¬ 
teresting to examine the effects of polydisperse crowd¬ 
ing agents on the complex problem of RNA folding to 
further some of our predictions. 


CONCLUSIONS 

In summary, using explicit simulations of crowding 
particles, we predict multiple and unexpected scenar¬ 
ios for the effects of polydisperse crowding environ¬ 
ment on the size and shape of a semiflexible polymer, 
which has served as a model for DNA and even RNA. 
Depending on the size, shape, and composition of the 
mixtures of crowding particles we find evidence for 
both compaction, and surprisingly dramatic increase 
in size as well. The results are of great relevance to the 
recent explosion of interest in the behavior of RNA [8], 
DNA [2BJ05], and proteins 149 > [50] in macromolecular 
crowding conditions both in vitro and in vivo. The 
prediction that shape of chains, such as DNA, RNA, 
and F-actin m, can be dramatically altered in a poly¬ 
disperse milieu can be tested in experiments. 


METHODS 

Model. To study the effects of crowding particles 
on a stiff chain, we used a coarse-grained model of 
WLC polymer (N m = 300), SS and SC crowding par¬ 
ticles (Figure la). The length scales, cr sp h , cr cy u and 
<r m (« 3.18 nm for DNA) denote the size (diameter) 
of the SS, SC crowding particles and the monomer of 
WLC polymer, respectively. We set the aspect ratio 
of the SC to be 2 and a sp h = 2 2 / 3 a cy i = 4cr m , so that 
the volumes of the individual SS and SC crowding 
particles are identical. In the WLC model, chain 
connectivity, with a fixed bond length, is maintained 
using a large spring constant connecting two consec¬ 
utive beads. The bending rigidity of the chain was 
implemented by quadratic bond angle potential. We 
chose Weeks-Chandler-Andersen (WCA) potential for 
interactions between monomers, and r -12 soft-sphere 
potential is employed for excluded volume interac¬ 
tions for crowder-crowder and crowder-polymer beads. 

Soft Spheres (SS) The energy function for the sys¬ 
tem consisting of the WLC and soft spherical crowders 
is, 


E = Es + Eb + F/wca + Er, 

(2) 

Es = N £k^ -«l-<•>’, 

i=1 L ° 

(3) 

N m - 1 


bd 

to 

II 

M 

1 

0 

to 

(4) 


2=1 
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FIG. 7: (a) A schematic of the expected changes in the size of biopolymers in the milieu of E. coli -like poly disperse 

crowding environment with cj) = 0.3 as L/l p is varied. Depending on the parameter value 7 = L/Z p , which characterizes 
the chain length and stiffness, the polymer undergoes swelling (7 > 10 ) or coil-to-globule transition (7 10 ). (b) 

PEG-induced compaction of Azoarcus ribozyme in 0.56 mM-Mg 2+ ion solution (blue circles). For comparison we show in 
simulation (black circles) results for R g (4>) changes in monodisperse SS crowding from^. 


Ewca = ^ © 

bj m 





^ = £d- 

v r « 


12 


a i + a a 


? u 13 


(6) 


where cr^ = , 0 (x) is a Heaviside step function, 

N m is the number of monomers, K (> 1000 &bT//q) 
is the spring constant with Z 0 being the bond length. 
The bending rigidity constant is G, cq is the diameter 
of a bead, 6i is the angle between the monomer bond 
vectors {fi +1 — g) and (rq — G-i), and e is Lennard- 
Jones energy constant controlling the strength of the 
excluded volume interaction. The diameter of the SS, 
<j sp h = 2 dm where cr m is the size of the monomer. We 
chose N m = 300 in our simulations. 


Crowders as spherocylinders (SC). We model 
the spherocylindrical crowder by connecting five 
spherical crowders allowing for overlap (Figure la). 
Five beads in the anisotropic crowders are connected 
using Eg and Eb in Eqs. © and © with very 
large values for the spring and the bending rigidity 
constants ( K sp and G sp , the analogues of K and G in 
[3] and [4| in order to maintain the cylindrical shape. 
We ignore excluded volume interaction between 


the beads within a particular cylindrical crowder, 
because the parameter G sp is sufficiently large. By 
choosing the diameter of the cylinder a cy i = 2.36cr m , 
the volumes of SS and SC crowders are identical 

(h a sph = h°lyi + x 2(7 cyl) (Figure la). The 

parameter values are given in Table SI. 


Mixture of SS and SC. To examine if the crowders 
of different shape have additive effects on the size 
of the WLC, we also considered a system containing 
spheres and spherocylinders. We chose equimolar 
mixture containing Ni = N 2 = N/2 spheres and 
spherocylinders so that the total volume fraction of 
the crowders is 0 -= N\ (vss + vsc ) = Nvss with 
vss = vsc • Thus, in a mixture fiss = 4>sc = 0/2 
where (j)ss and (j)sc are the fractions of volume 
occupied by SS and SC, respectively. 

Modeling the E. coli environment. Using the 
approximate composition of the crowding particles in 
E. coli in terms of sizes of crowders m we mimic 
the cytoplasm as a mixture of spheres containing 
the three largest particles. They are the ribosome 
with radius r\ = 10.4 nm, polymerases and other 
large complexes with average size r 2 = 5.2 nm, and 
smaller complexes where mean size r% ~ 2.6 nm. The 
composition of the three classes of particles is 11 %, 
11 %, and 8 % respectively. These simulations provide 
a general framework for understanding the fate of stiff 
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molecules in cell-like environment. 

Simulation details. In order to obtain adequate 
sampling of the conformational space of the system, 
we used low friction Langevin dynamics (LFLD). It 
can be shown rigorously, and has been confirmed in 
simulations, that the thermodynamic properties of the 
system do not depend on the choice of the friction coef¬ 
ficient [7 . In the LFLD, the diameter of the monomer 
<r m , r = (mcr^/e) 1 / 2 , and e were chosen as the units 
of length, time, and energy, respectively. The value 
of crm suitable for DNA is ~ 3.18 nm. Friction coef¬ 
ficients, £ m , for monomers and ( c for crowders were 
C m = 0.05 mr _1 and Cc = Cm^sph/^m [52] (Table SI). 
The duration of each trajectory ranges from 2 x 10 7 to 
5 x 10 7 At where At = 0.01 r. 

Initially, a semi-flexible chain was placed in a sim¬ 
ulation box without the crowders. We performed the 
LFLD for 10 7 time steps to equilibrate the system. 
The crowding particles were added to generate a sam¬ 
ple with cj) = 0.05. Higher volume fractions were 
reached by inserting additional crowding particles to 
the simulation box. Subsequently, Lennard-Jones in¬ 
teraction annealing (adiabatic increase of e) was car¬ 
ried out for (j) > 0.3 in order to improve the speed of 
equilibration, and to avoid catastrophic crashes during 
the insertion of particles. To be specific, at the begin¬ 
ning of annealing, we inserted the crowding particles 
at random positions by assigning e = 0.1 ksT and 
O.lAt for the integration time step of LFLD. First, the 
Lennard-Jones interaction parameter e was increased 
by 0.05 ksT for every 10 4 time steps until it reaches 
e = 1.67 &bT. Next, we increased the simulation time 
step by 0.1 At for every 10 4 time steps until the time 
step reaches At. 

For purposes of efficient computation we devised a 
method (see also SI in Ref. m) in which crowding 
particles are added on the fly and the volume of the 
simulation box is adjusted to keep the volume frac¬ 
tion constant. During the simulations we adjusted the 
size of the simulation box according to the chain con¬ 
formation at a given time to minimize the number of 
crowders. At every time step, we checked if the chain 
is enclosed in the simulation box. If any monomer and 
the boundary of the box is closer less than three times 
the average distance between the crowders, we resized 
the box and added crowding particles to the newly ex¬ 
tended empty spaces. As a result of constantly resizing 
the box (a cuboid with changing dimensions) the num¬ 
ber of crowders varies. The volume of the cuboid and 
the number of crowders are varied in such a way that 
0 is a constant. The average number of crowders in 
our simulations varies from 4000 to 8000 depending on 
(j). 


The particular ensemble used in these simulation is 
not used frequently although it is discussed by Callen 
(see page 148 in [53]). If the system consists of mono- 
disperse particles (generalization to multi-component 
system follows readily as explicitly shown in Ref[53]) 
in which the number of particles ( N ) and volume 
(V) fluctuate, the independent variables conjugates to 
these two variables are chemical potential (/i) and pres¬ 
sure (p). The thermodynamic potential in this ensem¬ 
ble is p,T) = U - TS + pV - fiN = G - pN = 0 
(follows from Euler relation), which means that vari¬ 
ations in N and V do not change the thermodynamic 
potential fl(/qp, T). This is precisely what is desired 
in these simulations. 

To ensure that the results do not depend on the 
choice of ensemble we also repeated the simulations 
in the canonical ensemble for the polydisperse case. 
As expected on theoretical grounds the results for R g , 
P(R g ), and the energy per particle (thermodynamic 
quantity) in the two ensembles are the same. The 
comparison is given in the Fig.S6. 

In total, we generated 25 trajectories at each 
volume fraction to obtain statistical properties. We 
collected data for analysis after a minimum of 10 6 
simulation time steps. 
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SUPPORTING INFORMATION 


Multiple layered neighbor list. The large sys¬ 
tem size needed to reliably simulate the WLC chain 
in the presence of explicit crowders is computationally 
demanding. To To circumvent the system size prob¬ 
lem, we develop a computer code that minimizes the 
number of operations to compute interaction poten¬ 
tial. In order to acheive this goal, we devised and im¬ 
plemented the multiple layered neighbor list (MLNL) 
technique. To our knowledge, this methodology has 
not been used in simulations before. In conventional 
Verlet algorithm, the list of neighbors, which are the 
particles located within a cut-off distance, R c , is cre¬ 
ated for each particle at the beginning of the simu¬ 
lation. When a pair-interaction is needed, we search 
only the neighbors instead of computing interactions 
between all pairs of particles. Thus, the computa¬ 
tional cost decreases as R c gets smaller. However, 
since the positions of the particles are constantly evolv¬ 
ing, we have to update the neighbor list with a cer¬ 
tain frequency. In conventional Verlet list, all neigh¬ 
bor lists are updated whenever the maximum displace¬ 
ment of any particle exceeds R c . The frequency of up¬ 
dates increases as R c becomes smaller, thus increas¬ 
ing the computational costs for updating the neighbor 
list. The two competing demands (frequent update for 
small R c and infrequent update for computations in¬ 
volving larger number of interaction pairs) requires an 
optimal value of R c . 

MLNL is designed to reduce the computational costs 
for calculating interaction potentials and updating 
neighbor lists by using multiple numbers of neighbor 
lists. It consists of several neighbor lists each with 
a different cut-off distances, R^ < R ^ ••• < R ^. 
Interaction potentials are only calculated using the 
upper-most layer, which has the smallest R c = R ^. 
This strategy minimizes the cost of calculating inter¬ 
action potentials. When the maximum displacement 
of a particle exceeds R^\ instead of calculating 
the distance between all pairs of particles as in the 
conventional algorithms, we update the upper-most 
neighbor list using the neighbor list with R K C ’. Thus, 
the requirement of computing O (V 2 ) interactions 
is avoided to a large extent by using this technique. 
The disadvantage of the MLNL is that memory 
requirement can be quite large especially when the 
system size is large. For N = 300, three layers suffice. 
With this choice we were able to perform converged 
simulations. 



FIG. SI: Asphericity (A) and shape parameter (5) of 
WLC as a function of SS and SC volume fraction, 0. 


shape anisotropy of polymer. The asphericity param¬ 
eter (A) and a shape parameter (S'). Both A and S, 
which are rotationally invariant, are defined using the 
inertia tensor, 


T a f3 = 


1 

2 JP 


N 

hj 


- r ja ) (r i0 - rjfj) 


(SI) 


where ri a is the a (= x, y, z) component of bead i of 
the WLC chain. The eigenvalues of T a p are related to 
R g via R 2 = Tr(T) = JV Xj . The anisophicity (A) 
and shape (S) parameters are 



i=1 


(A,-A ) 2 
(Tr(T)) 2 ’ 


27 nL (A» - A) 

(Tr(T)) 3 


(S2) 


where A = ■ Tr ^ T ^ . For a globule, A = S = 0. Thus, 
we expect that if crowders induce compaction then A 
and S should decrease monotonically as 0 increases. 

Compaction in R g for 0 < 0.3 can also be correlated 
with shape changes. We calculated asphericity (A) 
and shape parameter (S) as a function of crowder 
characteristics (Fig|Sl|). For spherical crowders, both 
A(0) and 5(0) decrease with increasing 0, implying 
that the WLC becomes more spherical with increasing 
compaction. In contrast, there is a sharp increase 
in A (cj)) and 5(0) as the volume fraction of the SC 
crowders increases beyond 0 ~ 0.2 (FigJSl]). The 
value of A(0.4) ~ A(0) and 5(0.4) ~ 5(0), which 
shows that at both 0 = 0 and high 0 the WLC is an 
anisotropic ellipsoid (5 > 0). The elongation of the 
chain along the local nematic axis is also reflected in 
an increase in 50(0) when 0 > 0.2. 


Shape anisotropy of polymer. Asphericity (A) 
and shape parameters (5) are used to characterize the 


Analytical expression for the end-to-end dis¬ 
tance distribution of WLC. The end-to-end dis¬ 
tance distance distribution of WLC, obtained in Ref. 
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C m 

cr c 

K 

G 

K-sp 

G sp 

lo 

do 

k B T 

At 

C m 

Cc 

3.18nm 

4 dm 

1500e 

4.91e 

3000e 

15e 

1.11 cr m 

0 

0.6 e 

0.01 r 

0.05 mr” 1 

Rs) 


TABLE SI: Parameters characterizing the model. Lennard Jones energy constant e, the diameter of monomer cr m and 

/ 777,0-2 

r = y are used as the fundamental units for energy, length and time scales. K and G define the strength of bond 
and angle potentials (Eqs. 2 and 3); K sp and G sp are the corresponding parameters for SC crowders. lo is a bond length 
between monomers of a chain, k B T is a temperature, At is a simulation time step, Cm and Cc are the friction coefficients 
for monomer and crowders, respectively. 



FIG. S2: a. The time-dependence of the intermediate 
dynamic scattering function (Eq|sg. b. Energy metric 
(Eq |S5| ) as a function of t. We calculated these quantities to 
ensure that F q (t) —> 0 and d(t)/d( 0) —> 0 (or d(0)/d(t) ~ t 
in the inset). Because the averages are computed using t > 
10 7 At, whereas measure of ergodicity decay on t « 10 4 At, 
the chain is fully equilibrated at = 0.3 in Figure 5. 


P(r) 


47tA/V 2 

(1 - r 2 ) 9/2 


exp 


3 t 
4(1— r 2 ) 


(S3) 


with r = Ree/L, t = L/l p , U = 4 ^ 2e 3 -i+i 5 c- 2 and 
c — \t where L = (7V m — 1 )a m is the contour length 
was used to obtain numerical values of the persistence 
length, l p . 


Pair correlation function of polydispersed 
crowders. Homogeneously distributed, the average 
separations (D) between crowders in cell lysate would 
be Dx/cr x = (47 t/3 ) 1 / 3 (1/2)0- 1 / 3 =1.68, 1.68, 1.87 


for X = 1, 2, and 3, respectively [15]. In contrast 
to this expectation, the pair correlation between the 
polydisperse crowders of each size indicates local size 
ordering, a major population at r/a x ~ 1, which 
is more pronounced than in the monodisperse case 
(Fig{S4|. 


ing particle sizes are large care must be taken to ensure 
that the system is well equilibrated. To provide ex¬ 
plicit evidence that our computational results (0 = 0.3 
case in Figure 3 in particular) are obtained from full 
equilibration, we evaluated multiparticle correlation 
function, which enables us to discern a glass-like dy¬ 
namically arrested state from an fluid-like equilibrated 
state. The intermediate dynamic scattering function, 
which probes the density-density correlation in Fourier 
space at q = | q\ = 2i:/r s where r s is the position of 
the peak in the total pair distribution function, tells us 
how our system relaxes from an initial configuration: 

1 N 

F ? (t) = — ^ (S4) 

V k=l 

where fk(t) is the position of k- th monomer of our 
semiflexible chain at time t. As shown in Fig|S2^i, 
Fq(t) 0 on time scales that are nearly three orders of 
magnitude less than the length of trajectories used to 
affect equilibrium data. Thus, the memory of the ini¬ 
tial configuration is fully erased in our simulations. In 
addition, as an alternative check for the equilibration, 
the energy metric for our system [54ff56j . 

1 N _ _ 

d(t) = -L(4,,(i)-%(t)) 2 (S5) 


where E a ^(t) = \ f* E a ^(r)dr is the energy of particle 
i averaged over time t from the trajectory generated 
from two different initial condition a and /3, shows 
that d(t)/d( 0) -T 0 rapidly (Fig{S2|p) . Furthermore, 
d(0)/d(t) ^ £, a hallmark of ergodicity, is explicitly 
shown (Fig|S2]3 inset). The results show the stiff chain 
in the polydisperse crowding environment ergodically 
explore the conformational space, thus ensuring that 
the results in Figs.5 and 6 represent converged results. 


Evidences of equilibration. Because the crowd- 
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FIG. S3: Distributions of gyration radius of SAW chain 
in the absence of crowders (0 = 0.0) and polydispersed 
crowder environment (0 = 0.3) mimicking the cytoplasmic 
condition of E. Coli. Note that unlike WLC in the main 
text (Figure 5), SAW chain is compacted in the polydis¬ 
persed crowding environment. 



FIG. S4: Crowder-crowder pair correlation functions be¬ 
tween crowders with the same kind in the polydisperse 
crowding solution. 
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FIG. S5: Effect of poly disperse crowding environment on 
the size of polymers with different 7 = L/l p 




FIG. S 6 : Comparison of simulation results of WLC under 
polydisperse mixture environment in /xPT and NVT en¬ 
sembles. a. ^-distributions of WLC in NVT (circles) and 
li pT ensembles (red line), b. Average energy per particle 
calculated for the two different ensemble. The value for 
each particle with different size is identical in the two en¬ 
sembles, indicating that simulation results do not depend 
on the choice of ensemble. 











